A symmetry breaking mechanism for epithelial cell polarization 
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In multicellular organisms, epithelial cells form layers separating compartments responsible for 
different physiological functions. At the early stage of epithelial layer formation, each cell of an 
aggregate defines an inner and an outer side by breaking the symmetry of its initial state, in a pro- 
cess known as epithelial polarization. By integrating recent biochemical and biophysical data with 
stochastic simulations of the relevant reaction-diffusion system we provide evidence that epithelial 
cell polarization is a chemical phase separation process induced by a local bistability in the signaling 
network at the level of the cell membrane. The early symmetry breaking event triggering phase sep- 
aration is induced by adhesion-dependent mechanical forces localized in the point of convergence of 
cell surfaces when a threshold number of confluent cells is reached. The generality of the emerging 
phase separation scenario is likely common to many processes of cell polarity formation. 



I. INTRODUCTION 

The development of epithelial tissues (e.g., kidney 
tubules, respiratory and gastrointestinal tracts, etc.) re- 
sults from complex morphogenetic processes implying 
the arrangement of cells in layers organized along spe- 
cific directional axes [TJ [2] • Epithelial cells are endowed 
with a self-polarization mechanism defining an 'inner' 
and 'outer' side, which is mandatory to allow organs to 
exert their vital functions. In a well established in vitro 
cell system, which recapitulates the in vivo morphogen- 
esis, after a single epithelial cell is seeded in a three- 
dimensional gel (Fig. [TJi) , cell division begins, and a mul- 
ticellular aggregate arises [TJ. The cells in the aggre- 
gate are bound each other (through cadherin molecules, 
Fig. [Tb) and to an extracellular matrix (through integrin 
molecules, Fig. [ft). When the cell number reaches 5-6 
units, an inner cavity, named lumen, is spontaneously 
opened [3] (Fig. [T^,b). Afterwards, cells develop a top 
(called apical) and a bottom (basolateral) side (Fig. [TJi) 
having different chemical features, while cell-cell and cell- 
matrix contacts only persist in the basolateral region. 
Finally, the border between apical and basolateral sides 
is sealed by ring-shaped tight junction proteins, which 
spontaneously find their functional position (Fig. [TJi) and 
prevent intermixing of chemical components between the 
apical and the basolateral membranes, as well as the out- 
pouring of liquids from the lumen. The full polar- 
ization process has a complex nature and involves dif- 
ferent factors and stages. However, recent experiments 
have determined its master regulator (see [H [5] and ref.s 
therein): intracellular asymmetry and lumen opening 
are controlled by PIP2/3 phospholipids (see Fig[2| and 
their interaction with PTEN/PI3K enzymes which in- 
duce PIP2/3 segregation to opposite poles, while the 
PAR complex further stabilizes axial polarity (see the 
review |Tj and ref.s below). Even in such an in vitro sys- 
tem, however, the mechanism whereby the cell original 
spatial symmetry is spontaneously broken and polariza- 



tion develops remains mysterious [TJ- 

To describe cell differentiation, polarization and signal 
localization mechanisms, stochastic reaction-diffusion [5] 
and coupled kinetic rate equations have been widely used 
for intracellular signaling, gene regulation and autocat- 
alytic reaction systems (see e.g. [THS] and ref.s therein). 
In this context, here we use reaction-diffusion equations 
to model the PIP2/3 master regulator in order to ad- 
dress some open questions £Q, namely the mechanisms 
for: i) the lumen site choice; ii) its opening; Hi) the 
control of its final size; iv) the localization of tight junc- 
tions. After delineating the chemical reactions involved 
in the process, by simulation of its master equation we 
show that self-polarization can be understood in terms 
of statistical physics concepts as a symmetry breaking 
mechanism driven by the chemical regulatory network. 
We finally interpret the simulation results in terms of a 
simple mean-field model. 



II. MODEL 

Recent experiments have shown that PIP2 and PIP3 
for a module that acts as a master regulator control- 
ling all signaling pathways and cytoskeletal dynamics re- 
quired for epithelial cell polarization [H [5] . PIP2 /PIP3 
levels are regulated by the counteracting enzymes PI3K 
and PTEN, which respectively catalyze the switch of 
PIP2 to PIP3 and vice versa [TJ [10] (Fig. [2|>. The phos- 
pholipids (PIP2/3) are stably localized in the inner face 
of the cell membrane where they diffuse. The enzymes 
(PI3K/PTEN) diffuse instead in the cell volume, where 
they are present in limiting amounts, and become ac- 
tive upon association with membrane spanning proteins 
or lipids. PTEN levels in the membrane are controlled 
by its binding to PIP2, thus realizing a positive feed- 
back loop (see Fig. |2J). PI3K levels in the membrane 
are controlled by its binding to cell-cell adhesive recep- 
tors (cadherins) and cell matrix adhesive receptors (inte- 
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FIG. 1: (Color online) a,b) After a few cell divisions, at the 
common point of contact among the cells of an aggregate the 
intrinsic high membrane curvature locally induces formation 
of a microlumen (inner circle) and the breaking of cell-cell ad- 
hesive receptor links (pink and purple bars). This breaking fa- 
vors the appearance of a small region rich in PIP2 lipids (blue) 
which provides a nucleation center for the ensuing symmetry 
breaking process. b,c) Everywhere but in correspondence of 
the PIP2-rich region, the cell is surrounded by adhesive con- 
tacts, either cell-cell (cadherins) or cell-matrix contacts (inte- 
grins). Cad (pink bars) and Cad* (purple bars) respectively 
indicate inactive and active cadherins. d) Scheme of the api- 
cal (blue) and basolateral (red) regions and the ring of tight 
junctions (green) spontaneously assembling at their interface. 
The extracellular matrix (yellow) is indicated as well. 



grins, schematically indicated by C/M in Fig. [2]) [TT] . 
To bind PI3K, cadherins must be activated (Cad* in 
Fig. [TJa and [2| by engagement with cadherins of a neigh- 
boring cell (named C/C in the diagram of Fig. [2]) PI3K 
is active only when associated to either activated cad- 
herins or intcgrins. Since PIP3 stabilizes the activated 
form Cad* [T^] , these interactions create a positive PIP3- 
PI3K feedback loop, mediated by the existence of cell- 
cell contacts (Fig. [2]). Before polarization, cadherins and 
integrins are activated along the whole membrane and 
PIP3 uniformly prevails on PIP2 determining a stable 
PIP3-rich phase over the whole membrane. A local de- 
pletion of PIP3-PI3K can be created if a large enough 
membrane area with disrupted cell-cell links is formed, 
thereby breaking the PIP3-PI3K feedback loop (Fig. [2]) 
and originating a germ of a PIP2-rich phase (Fig. lip 
and[2|. Then, the PIP2-PTEN feedback loop may locally 
prevail, inducing a PIP2 and PIP3 surface compartmen- 
talization that splits the cell membrane in two regions, 
or phases [T3J, characterized by different chemical con- 
centrations of the signaling molecules. 

Several mechanisms have been proposed j3] to explain 
the initiation of the polarization process, which is typ- 
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FIG. 2: Reactions scheme of the epithelial polarization net- 
work. The core players are the PI3K and PTEN enzymes, 
which respectively catalyze the switch of the PIP2 lipid into 
PIP3 and vice versa. PTEN becomes effective upon associa- 
tion to PIP2, while PI3K becomes effective upon association 
with activated cell-cell adhesive receptors (cadherins, Cad*) 
or extracellular matrix adhesive receptors (C/M). To bind 
PI3K, cadherins (Cad) must be activated (Cad*) by linking 
other cadherins on neighboring cells (C/C) and by the ac- 
tion of PIP3. The symmetry of PTEN/PI3K in generating 
PIP2/3 can be broken by cadherins inactivation. The PIP2- 
PTEN feedback loop may then locally prevail, originating the 
chemical polarization of the cell. 



ically found in high curvature membrane regions, es- 
pecially in correspondence of multiple cell-cell contacts 
(cellular vertex) [2] away from the extracellular matrix 
(Fig. [I]). We observe here that cell-cell contacts can be 
broken [T3] as soon as the adhesion energy W a (i.e. the 
energy per unit area needed to break cadherin contacts) 
becomes comparable to the elastic energy stored in the 
cell membrane in high curvature regions, such as at the 
confluence of several cells (Fig. [T^,,b). Since the elastic 
energy per unit area W c stored in the membrane is |15f - 



W c = n/2r 



(1) 



where k is the membrane bending rigidity and r the local 
curvature radius [IB], the condition W a ~ W e allows to 
estimate the critical curvature radius r a where cell-cell 
contacts start being disrupted as [T7] 



(2) 



The critical value r a can be easily estimated. In eu- 
karyotic cells, the typical adhesion energy of cadherin 
contacts is W a ~ 10 _11 /iJ//^m 2 |17j . while the typical 
bending rigidity is k ~ 400fc B T ~ 16 x KT 13 iiJ Q7], 
giving [37] 



0.3 ^m . 



(3) 



According to the above estimate, when the number of 
cells in the initial aggregate increases up to the 5-6 cell 
stage, at the cell convergence points (see Fig. [Tp) the 
membrane curvature increases as well and, especially in 
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areas not in contact with the extracellular matrix, cad- 
herin bridges are subject to forces that can disrupt links. 
By such a mechanism, a local opening of cell-cell contacts 
breaks the PIP3-PI3K feedback loop and induces a local 
unbalance towards PIP2 formation (Fig. [2]) and a germ 
of the PIP2-rich phase can be ushered in (Fig. [T]b) . 

In Sect. |III| we show that only germs of the PIP2-rich 
phase larger than a treshold radius r t hr actually survive 
and grow. This fact suggests that although the uniform 
PIP3-rich phase is not the more stable state for the sig- 
naling network, a polarized state characterized by the co- 
existence of the PIP3-rich and the PIP2-rich phase may 
be reached only by overcoming a barrier in a suitably 
defined effective energy (see Sect. IV for a detailed dis- 
cussion of this point). Therefore, we are faced with the 
following physical picture: if elastic forces due to high 
membrane curvature in the region of cell convergence 
(cellular vertex) trigger disruption of cadherin links in a 
region of size larger than fthrj the PIP2-rich patch grows 
favoring further breaking of cadherin links (Fig. [2]) and 
the formation of a lumen. Thus, for the process of lumen 
formation to start, it is necessary that the local curvature 
radius in the cellular vertex satisfies 



r a > r t hr 



(4) 



The growth of the PIP2-patch and lumen slows down 
and eventually comes to a stop as soon as cytosolic 
PTEN is depleted. This way, at the end of the pro- 
cess the cell reaches a stable polarized state characterized 
by the coexistence of the PIP2-rich and the PIP3-rich 
phase [LSI [19] . and a lumen coinciding with the PIP2- 
rich phase is formed. 



III. SIMULATIONS. 

In this Section we investigate on quantitative grounds 
the above described scenario of polarization. 

Since the chemical reaction and diffusion processes are 
intrisically noisy, we simulate the corresponding dynam- 
ics by a stochastic algorithm, using realistic values for 
reaction and diffusion rates. We can check this way that 
noise alone is not sufficient here to overcome the energy 
barrier separating the uniform and polarized state in ob- 
servational times, if an initial PIP2 seed of size larger 
than rthr is not created by an external interaction. 

We represent the plasmamembrane by a lattice of 
TV = 10242 (mostly hexagonal) sites of area a ~ S/Af ~ 
(0.1 /jm) 2 on a sphere surface S = AirR 2 with radius 
R = 5 /im (Fig. [3]) . Each site is populated by a number 
of molecules of the chemical factors and their dynamics 
is described by standard master equations. For instance, 
the PIP2 — > PIP3 process is described by 

d t P(N PIP2 ,N P1P3 ,...) = 
+ VF(PIP 3 ->• PIP 2 )P(^pip 2 - 1, iVpip 3 + 1, . . .) 
- VF(PIP 2 ^PIP 3 )P(/Vpip 2 ,/Ypip 3 ,...), (5) 
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^6^PI3K* 


k 6 = 0.1 


Cad -> Cad* 
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k 7 = 2 ■ 10~ 5 


Cad* — > Cad 




fc 8 = 0.5 



TABLE I: List of chemical reactions involved in epithelial 
polarization and their corresponding rates. X* denotes the 
membrane-bound, activated form of molecule X. Rate con- 
stants k are given in s _1 . Activation rates can be trans- 
formed to (sM)" 1 units multiplying by Avogadro's number 
Na- Michaelis-Menten constants K are pure numbers and 
can be transformed to M units dividing by Na/N- For the 
rate constants values see references in |2U] . 



where P(N%, ■••) is the probability to have at time t a 
number Nx of type X molecules at a given site (say, i.) 
The list of relevant reactions with their corresponding 
rates W is given in Table [jj [35]. PIP2/3 diffusion is de- 
scribed by random jumps of a molecule from site i to its 
neighboring site j with rate W(i — > j) = NxD/a, where 
D = 0.5 /im 2 /s is phospholipid diffusivity [20]. Since the 
diffusivity of cytosolic enzymes (PI3K/PTEN) is much 
larger than that of membrane pospholipids [5D] , their dis- 
tribution in the cytosol is treated as uniform. 

For the simulations we use a variation of Gillespie algo- 
rithm [21J taking into account the spatial non- uniformity 
of the system. At time zero, a random number is gener- 
ated to determine the next reaction or elementary diffu- 
sion process to occur, with a probability proportional to 
the corresponding W factor from Tabic [T] Then, time is 
advanced as a Poisson process with rate again determined 
by the W factors. These steps are repeated iteratively 
until the desired simulation time is reached. 
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FIG. 3: (Color online) Growth of the PIP2-rich phase (blue 
upper patch). The color scale shows the gradation of PIP 2 
content: ip is the relative concentration difference between 
PIP3 and PIP2 at a given site. The system at initial time 
is in a uniform PIP3-rich phase, apart from an initial PIP2- 
rich seed germ of size ro larger than the threshold radius r t hi 
(small circle.) After 4 minutes a PIP2-rich patch becomes ap- 
parent and its radius saturates after approximately 10 min- 
utes to the equilibrium value r eq . 
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FIG. 4: (Color online) Growth of the PIP2-rich phase induced 
by a germ of initial radius ro. Germs with ro smaller than a 
threshold radius r t hr ~ 0.3 fim are melted by diffusion, while 
larger germs grow to the equilibrium value r cq . Error bars 
show standard deviations computed over n — 50 different 
random realizations of the process. 



We suppose that a circular PIP2-rich patch of radius ro 
is initially formed in the sea of the PIP3-rich phase [39] 
and investigate its dynamics to check whether a stable 
polarization state is attained (Fig. [3|. Fig. [4] shows the 
time evolution of circular patches of different initial radii 
ro- Patches smaller than a threshold radius rthr ~ 0.3 ^m 
are dissolved by diffusion and thermal processes, and do 
not impair the stability of the uniform PIP3-rich phase. 
Conversely, patches larger than r t hr grow in time trigger- 
ing the separation of the cell surface in a PIP2-rich and 
a PIP3-rich region and eventually reach an equilibrium 
(Fig. [3]). Notably, the threshold radius r t h r ~ 0.3 /zm 
derived from the above calculation is consistent with the 
previously independently derived value for the adhesion 
radius, r a . The two phases are divided by an interface 
of characteristic width 5 ~ \J D jk c ~ 1 /Ltm [T^] , where 
k c is of the order of the catalytic constants of the two 
catalytic reactions of Table [TJ (first two rows). 

The kinetic of this heterogeneous nucleation process 
can be understood in terms of non-equilibrium, reaction- 
diffusion stochastic dynamics. In reaction-diffusion sys- 
tems instabilities are often produced by Turing's mecha- 
nism [22j . Here we find that pattern formation starting 
from a locally stable homogeneous state is triggered by a 
local perturbation by a nucleation center of size ro larger 
than a critical size r c [T3] . 

Fig. [5] shows that the equilibrium size r eq of the PIP2- 
rich patch, and therefore of the lumen, is controlled by 
the number of PI3K and PTEN molecules. In the absence 
of any limiting mechanism, the growth of the PIP2-rich 
patch would in fact lead to a PIP2-rich phase completely 
invading the cell surface. However, due to the coupling 
to a finite PTEN and PI3K reservoir, the system self- 
tunes to a phase-coexistence state and the process stops 
when the PIP2-rich patch reaches the equilibrium size 
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FIG. 5: (Color online) The equilibrium lumen size r oq is 
an increasing function of the relative number of PI3K and 
PTEN molecules. Inset: radial distribution of PI3K (red 
squares, jVpisk = 1.26 ■ 10 s molecules) and PTEN (blue cir- 
cles, ./Vpten = 0.43- 10 5 molecules), from the patch center, at 
equilibrium. 

r eq 18 20 . Interestingly, the fact that the size of the 
PIP2-rich patch, and consequently of the lumen, is con- 
trolled by the precise number of PTEN molecules, is in 
qualitative agreement with the observation that deletion 
of a single PTEN allele can interfere with the polarization 
process [25] . 

The observation in the present scenario of a compar- 
atively large threshold radius, of the order of one tenth 
of the cell size, suggests the existence of a correspond- 
ingly large barrier of effective energy dividing the uni- 
form state from the phase-separated one. This prevents 
thermal and chemical noise from triggering spontaneous 
symmetry breaking and lumen formation. However, an 
external mechanical action creating a sizeable PIP2-rich 
patch, due to the presence of localized regions of high 
membrane-curvature, can overcome the barrier and start 
polarization [30]. Our picture also explains tight junc- 
tion localization. Experimental data show that the sta- 
ble binding of tight junction proteins to the membrane re- 
quires both a protein complex named PAR3-PAR6, which 
is localized in the PIP2-rich phase by a chain of reactions, 
and cell-cell contacts, which are maintained only in the 
PIP3-rich phase [27]. The spontaneous aggregation of 
tight junctions is thus constrained by a biochemical logi- 
cal AND to take place only on the ring-shaped boundary 
separating the PIP2-rich from the PIP3-rich phase. 

IV. MEAN-FIELD 

In this Section we show that the results of the simula- 
tions can be conveniently interpreted in terms of an ef- 
fective mean-field model, following the approach detailed 
in Refs. QS[5S]. 

Fast diffusion of PI3K and PTEN enzymes in the cy- 
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tosol, and the conservation law [PIP2] + [PIP3] = c allow 
to effectively describe the state of the cell membrane in 
terms of the configuration of the single-component con- 
centration field [raiES] : 



cp = [PIP3] - [PIP2]. 



(6) 



The resulting effective equation for ip can be set in the 
simple Landau-Ginzburg form: 



d t <p = DV 2 <p + V \<p)+t, 



complemented by an integral costraint expressing the 
coupling of the concentration field ip to the reservoir of 
free cytsolic enzymes (see Ref.s [THl US] and Supplemen- 
tary Information in |20j.) In Eq. Q, D is the diffusivity 
of lipids on the cell membrane, V[(p) is an effective po- 
tential, and £ is a stochastic term taking into account the 
effect of thermal and chemical noise. 

The mean-field effective potential V{<p) can be easily 
derived, via a quadratic approximation, from the stochas- 
tic model described in Sect. |III| under the assumption 
that the cytosolic PI3K, PTEN and Cad fields are in 
approximate equilibrium with the membrane PIP2 and 
PIP3 fields, and therefore "slaved" to the ip field pll25] : 



V'(<p) = 



c 2 ~ V 2 



2K + c + ip 2K + c - (p 



where 



&2 h 
ki fc 5 



[pten; 



free ; 



k- 



— [PI3K] frco • [Cad] free . 
k 6 k 8 



, (8) 

(9) 
(10) 



The terms in the r.h.s. of Eq. ^ describe respectively 
conversion of PIP3 into PIP2 due to the action of PTEN, 
and conversion of PIP2 into PIP3 due to the action of 
PI3K activated by cadherins (Fig. [2]). The quadratic 
terms oc c 2 — (p 2 encode respectively the PIP2 — >PTEN 
and the PIP3->Cad^PI3K feedback loops (Fig. [2]). In 
particular, a' = when cadherin links are broken. 

In a wide region of parameter space around the real- 
istic parameter values from Table [IJ the effective poten- 
tial V(tp) is bistable (Fig. |6j for a detailed description 
of the bistability region, see Ref. [IS]). The two poten- 
tial wells in Fig. [6] correspond to a stable PIP2-rich and 
a metastable PIP3-rich phase, separated by an energy 
barrier AV. 

The mean-field model ( |7|8| and the bistability of the 
effective potential V(<p) provide an interpretation to the 
simulation results, showing that the stable polarized state 
characterized by the coexistence of PIP2 and PIP3 in 
complementary regions is separated from the metastable 
PIP3-rich phase by an effective energy barrier AV. Ac- 
cording to the theory of Landau-Ginzburg equation, 
PIP2-rich seeds larger than a critical value are bound to 
expand in the PIP3-rich sea with a velocity proportional 
to 1/AV [23]. 
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FIG. 6: (Color online) Graph of the effective potential V(jp) 
defined by l]8|9|10 l, where the cytosolic values [PTEN] frec ~ 
1.1 • 10 4 , [PI3KjLo ^ 3.7 ■ 10 4 , [Cad] frco ~ 5.4 ■ 10 3 
(molecules/cell) were computed from simulations of the 
stochastic model described in Sect. |IH| at time t = 30 s, after 
the fast equilibration of the corresponding cytosolic pools of 
enzymatic factors with the initial configuration of the ip field. 
The effective potential V(ip) has two minima, the left hand 
one corresponding to a stable PIP2-rich and the right hand 
one corresponding to a metastable PIP3-rich phase. The two 
phases are separated by an effective energy barrier. 



The cytosolic concentrations [PI3K] frcc , [PTEN] free , 
and [Cad] froc appearing in (JsJ) may be expressed as in- 
tegrals of the concentration field tp [55J. The re- 
sulting global coupling has the effect of driving dynam- 
ically the cell membrane towards an equilibrium polar- 
ized state where the PIP2-rich and PIP3-rich phases co- 
exist: the growth of the PIP2-rich phase "eats up" free 
PTEN molecules from the cytosol, decreasing AV until 
phase coexistence is reached [35J . This process may 
be understood via a simple physical analogy with the 
non-equilibrium process taking place during the liquid- 
vapor transition in a sealed vessel: there, the rise of the 
vapor pressure (which in our analogy corresponds to the 
number of cytosolic enzymes) provides a negative feed- 
back, slowing down the growth of the vapor phase and 
eventually leading the system to a state of phase coex- 
istence. The main difference between the two systems 
is that in the liquid-vapor transition a local conservation 
law holds for the particle field, while in the growth of sig- 
naling domains on the cell membrane the tp field satisfies 
only an approximate global constraint [25] encoded in the 
integral expressions for the coefficients a, a' appearing in 
~|9p| - 



It is worth observing here that while spontaneous po- 
larization in eukaryotic directional sensing [20] can be de- 
scribed in terms of an homogeneous nucleation process, 
whereby seeds of a PIP3-rich phase are created by ther- 
mal and chemical noise in the sea of the PIP2-rich phase 
and grow by a coarsening process [191 125j . the present 
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picture of epithelial polarization reminds instead hetero- 
geneous nucleation, i.e. a situation where the effective 
potential barrier AV is so high that spontaneous nu- 
cleation does not occur in typical observational times, 
and needs to be triggered by the introduction of a large 
enough nucleation germ. 

V. DISCUSSION 

We have shown that a simple symmetry breaking 
mechanism, informed with the recently discovered bio- 
chemical and biophysical details of the system, accounts 
for a wealth of morphogenetic processes in epithelial po- 
larization. The model makes specific predictions on the 
dependence of the threshold radius, lumen size, tight 
junction positioning and width, on the biochemical sys- 
tem parameters. Our results shed light on the role of 
PTEN as a tumor suppressor protein, whose expression 
levels are known to be critical to prevent the onset of 
cancer. In particular, our models predicts that by de- 
creasing the number of PTEN molecules, the lumen size 
should decrease, and for very low PTEN levels no lu- 
men at all should form. The experimental validation of 
the model could be performed by genetic manipulation of 
the amounts or activity of cadherins, PI3K and PTEN. 

We have also shown that curvature induced forces are 
a very plausible candidate for triggering the symmetry 
breaking process at the right time. This could be veri- 
fied experimentally by trying to induce localized forma- 
tion of a growing PIP2 patch and lumen by mechani- 
cally breaking adhesion bonds in localized regions of the 
membrane of epithelial cells surrounded by extracellu- 
lar matrix. Under these conditions, our model predicts 
that only breaking adhesion bonds in regions larger than 



rthr should induce the formation of a growing patch of 
the PIP2-rich phase, while smaller PIP2-rich patches, in- 
duced by breaking adhesive bonds on smaller regions, 
should shrink spontaneously 

Interestingly, the bistable PI3K-PTEN module here 
described plays also a key role in chemotaxis, where PI3K 
is initially activated by chemotactic receptors (see the 
review |29j and ref.s therein) rather than by adhesive 
receptors. While experimental evidences and our results 
suggest that epithelial polarization is induced by a nucle- 
ation center of the PIP2-rich phase generated by mechan- 
ical forces, the polarization of migrating cells is likely to 
be triggered by spontaneous fluctuations in PIP3/PI3K 
levels US]. 

PIP3 localization also regulates chemotactic polar- 
ization [19j [20] and cell spindle orientation 30 . The 
similarity underlying the mechanisms in these very 
different aspects of cell life hints to the possibility 
that phase separation phenomena might have a general 
role in the cell [HI [20] and in its nucleus [3lJ 132] , 
The principles emerging here could explain in a uni- 
versal way the deep analogies observed in a variety of 
cellular processes involving spatial polarity formation [1] . 
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